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ABSTRACT 



Context. Flare Hard X-Rays (HXRs) from non-thermal electrons are commonly treated as solely bremsstrahlung (free-free = f-f), 
recombination (free-bound = f-b) being neglected. This assumption is shown to be substantially in error, especially in hot sources, 
mainly due to recombination onto Fe ions. 

Aims. We analyse the effects on HXR spectra y(e) and electron diagnostics by including non-thermal recombination onto heavy 
elements in our model. 

Methods. Using Kramers hydrogenic cross sections with effective Z = Z^ff, we calculate f-f and f-b spectra for power-law electron 
spectra within both thin and thick target limits and for Maxwellians with summation over all important ions. 

Results. We find that non-thermal electron recombination, especially onto Fe, must, in general, be included with f-f for reliable 
spectral interpretation, when the HXR source is hot, such as occulted loops containing high ions of Fe (f-b cross-section oc Z*). The 
f-b contribution is greatest when the electron spectral index S is large and any low energy cut-off is small, because the electron flux 
spectrum F(E) emitting f-b photon energy e is oc (£ = 6 - VzT^ (Vz is the ionisation potential) and not oc (E = 6)"*^' as for f-f. The 
f-b spectra recombination edges mean a cut-off E^ in F(E) appears as an HXR feature at e = E^ + Vz, offering an E^ diagnostic. For 
thick target sources, the presence of E^ appears as edges in J'ie), not in /(e), but it is still detectable. Including f-b lowers the E(E) 
needed for prescribed HXR fluxes greatly in some cases; and even when small, it seriously distorts F(E) as inferred by inversion or 
forward fitting of J{€) based on f-f alone. 

Conclusions. The f-b recombination from non-thermal electrons can be an important contributor to HXR spectra, so it should be 
included in spectral analyses, especially for hot sources. Accurate results will require use of better cross sections than ours and 
consideration of source ionisation structure. 

Key words. Atomic processes; Sun:corona-Sun:flares-Sun:X-rays, gamma rays 



1. Introduction 

Ever since their first detection (Arnoldy, Kane and Winckler 
1968, Kane and Andersen 1970), flare hard X-ray (HXR) bursts 
(photon energies e > 10 keV or so) have been recognised as 
an important diagnostic of electron acceleration and propagation 
(e.g. Brown 1971, Lin and Schwartz 1987, Johns and Lin 1992). 
The large electron flux and power imply they play a substantial 
role in flare energy budgets and pose challenges for electron ac- 
celeration mechanisms (see recent reviews by, e.g. Vilmer et al. 

2003, Brown 2005, MacKinnon 2006). Recent copious high res- 
olution HXR spectral data from the RHESSI mission (Lin et al 
2002) have created the possibility of detailed reconstruction of 
source electron spectra (following Brown 1971) offering impor- 
tant constraints on the electron energy budget and acceleration 
processes (Plana et al. 2003, Conway et al. 2003, Massone et al. 

2004, Kontar et al. 2004, 2005, Brown et al. 2006). 

In inferring electron flux spectra F{E), the HXR radia- 
tion mechanism has always been taken to be f-f collisional 
bremsstrahlung of fast electron impacts with atoms and ions, gy- 
rosynchrotron and inverse Compton radiation being negligible at 
these energies for solar magnetic and radiation fields (Korchak 
1971). Though included for thermal electrons in hot (a few keV) 
plasma, f-b recombination radiation from non-thermals seems 
to have been assumed negligible other than in a preliminary 



Send offprint requests 
j ohn@astro .gla.ac.uk 



I.e. 



Brown, 



e-mail: 



study by Landini, Monsignori and Pallavicini (1973). In view of 
the importance of details in the photon spectrum 7(e) (photons 
sec^'kev"') for accurate reconstruction of F(E), we re-examine 
this assumption, and conclude (cf Mallik and Brown 2007) that it 
is not valid under some conditions, which quite commonly exist 
in some flare HXR source regions. 

It is not the intention of this paper to analyse precisely the 
theoretical recombination radiation spectrum from fast electrons 
under conditions (e.g. ionisation structure) for specific flares 
which are typically both inhomogeneous and time dependent. 
Rather we give approximate theoretical estimates of how im- 
portant it may be relative to bremsstrahlung under various lim- 
iting conditions. Specifically, we compare the two in the sim- 
plest, Kramers, cross-section approximations, for limiting cases 
of plasma ionisation. The recombination emission rate per elec- 
tron is very sensitive to the ionic charge, being oc Z'^Az (Kramers 
1923) per plasma proton for hydrogenic ions of charge Ze and 
number abundance Az. Thus the emitted f-b flux and spectrum 
depend strongly on the ionisation state, hence the temperature, of 
the plasma where the fast electrons recombine. In practice this 
will involve several ionisation stages of several target plasma 
species (since Z*Az may be large even for small abundance Az), 
which will vary along the paths of the electrons and be time de- 
pendent. 

The paper is organised as follows. In Section 2 we briefly dis- 
cuss relevant processes and the cross-section approximations we 
use, and obtain expressions for the total continuum photon spec- 
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tral contributions j(e) expected from an electron flux spectrum 
F{E) from f-f and from f-b processes. In Section 3, we com- 
pare these for a power-law F(E) with low cut off at £ < Ec and 
for a shifted power-law, and discuss implications for flare elec- 
tron spectra and energy budgets under several limiting plasma 
ionisation assumptions. In Section 4 we look at thermal and 
non-thermal components to show how the relative importance 
of each contribution depends on conditions in the flare by vary- 
ing parameters around those for a specific real event. Section 5 
discusses the efi'ect of including the f-b contribution on inverse 
problem inference of F{E) from y(e) while Section 6 summarises 
our conclusions and suggests directions for future work. Details 
of some of the equations are given in Appendix A. In Appendix 
B we discuss the total emission spectra from extended volumes 
for thin target, colhsional thick target and thermal cases. 

2. Free-free and free-bound emissivity spectra 

2.1. General considerations 



In this section, we discuss only local emissivities j(e) (photons 
cm"^ sec ' per unit e ). Relativistic and directivity eff'ects are 
disregarded (E, e «i nigC^) since the f-b/f-f ratio is largest at 
low E. Then, if target atom/ion type t has density n, and the fast 
electron flux spectrum is F(E) (electrons sec"' cm"^ per unit E), 
jie) for a colhsional radiation process is 



j(e) = Sj-,(e) = ^,n, f F{E)^(e, E)dE, 



(1) 



where dQ,/d£(€, E) is the relevant cross-section per unit e for 
target species t and the integral is over the range of electron en- 
ergies relevant to species t. 

2.2. Bremsstrahlung 

In the case of f-f (bremsstrahlung), dQ,lde(e,E) is essentially 
the same for any state of ionisation of an atomic species Z (Koch 
and Motz 1959), and the t summation in (1) need only be car- 
ried out over elements Z to give, for element abundances Az (by 
number relative to hydrogen), and total proton (pH-H) density Up, 



jeie) = Hp-LzAz ^ F(£)^(e, E)dE, 



(2) 



where dQszldeie, E) is the bremsstrahlung cross-section for el- 
ement Z and £■„;,„ = e since any free-free transition can only 
yield a maximum e - E. The bremsstrahlung cross-section per 
nucleus Z scales as and can be written 



dQE 



de 



~ eE~ 



q{e, E) , e<E 



(3) 



(and zero for e > E). Here a - e^/ftc is the fine structure con- 
stant and re = e^/nigC^ the classical electron radius, while q(e, E) 
is the ratio of the actual cross section to the Kramers cross sec- 
tion (Kramers 1923), which is the factor in front of q. While 
this is only a first approximation, not suitable for accurate abso- 
lute spectral inversion/reconstruction algorithms (Brown 2005), 
it will be adequate for the present purpose of comparing f-f with 
f-b emission, which we also treat in the Kramer's approximation. 
Then (2) and (3) give, for bremsstrahlung, 



jBie) 



r 



F(E) 



dE, 



(4) 



where 

= ^z(bz = ^zAzZ^ 



(5) 



is the heavy element correction for bremsstrahlung, with ~ 
1.6 for the solar coronal abundances we use - see later. 



2.3. Recombination Radiation 

The situation here is more compUcated. Firstly, 2-body radiative 
recombination (we neglect 3 -body recombination) of a free elec- 
tron of energy £ to a bound level m of energy -V(Z, i, m) in ionic 
stage / yields a photon energy e, which, apart from quantum un- 
certainty, is unique, namely: 



€ = E + V(Z,i,m). 



(6) 



That is, when a fast electron does recombine, all of its kinetic 
energy E plus V goes into a photon of that energy, in contrast to 
bremsstrahlung where photons of all energies e < £ are emitted. 

Furthermore, for each element Z, there is a range of Z-i- 1 dis- 
tinct ion stages i each with its own distinct set of energy levels 
(m) and a set of Z, /, ;7i-dependent recombination cross-sections. 
Thus recombination collisions of a mono-energetic beam with 
a multi-species plasma gives rise to a set of delta-function-Uke 
spectral features at all energies (6) corresponding to elements 
Z, ionic stages i and levels m . For a continuous electron spec- 
trum, this yields a continuum photon spectrum that is a sum of 
an infinite series of energy- shifted electron flux contributions. In 
contrast to bremsstrahlung it does not involve an integral over a 
continuum of electron energies. 

For a general plasma the basic particle type "f" onto which 
recombination occurs is level m of ion stage ; of element Z with 
recombination cross-section diff'erential in e for that t: 



dQ. 



'Rt 



de 



(e) = QRAE-e+V,), 



(7) 



where Qnt is the total radiative recombination cross-section for 
species t and 6(E') is the delta-function in energy such that 
6{E')dE' - 1 . Then the total recombination emission spec- 
trum for electron flux spectrum F{E) is 

Me) = UpLtA, ,j G«,(e, E)6{E - e + V,)F{E)dE 

i:,A,npQR,(e, e - V,)F(e - V,), (8) 

where At is the numerical abundance of species t relative to rip. 
The forms for Qr,, for general f, are complicated and have to be 
calculated numerically, as do the values of A, when individual 
ionisation states are considered. However, in the Kramers ap- 
proximation (with unit Gaunt factors) there is an analytic expres- 
sion for hydrogenic ions, which we will use to estimate djR/de 
compared with djs/de, namely, for recombination onto level m 
of the hydrogenic ion of element Z (Kramers 1923, Andersen et 
al. 1992, Hahn 1997) 



Qr = 



327r '74,.2 



3 V3a 'm^eE' 



(9) 



where X = nige'^/ltp- is the hydrogen ionisation potential. 

For an element in its highest purely hydrogenic ion state the 
emissivity spectrum would then be 



. 32n r^^^^ IFie-Z^x/m^) 
jRzye) = — — 2.„ 



3V3ar 



e - ZP-x/mP- 



(10) 
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Table 1. Elements with their coronal abundances and ionisation potentials at T » 10^ K 



Element 


Z 


A- 




A.Z" 


V, = Z^x (keV) 


H 


1 


1 


1 


1 


0.0136 


He 


2 


0.096 


0.384 


1.536 


0.0544 


C 


6 


3.57 X HH 


0.013 


0.463 


0.490 


O 


8 


8.57 X 10"" 


0.055 


3.511 


0.870 


Ne 


10 


1.07 X 10-" 


0.011 


1.071 


1.360 


Mg 


12 


1.33 X IQ-" 


0.019 


2.755 


1.958 


Si 


14 


1.27 X IQ-" 


0.025 


4.871 


2.666 


S 


16 


1.61 X IQ-' 


0.0041 


1.053 


3.482 


Ca 


20 


8.50 X 10"^ 


0.0034 


1.360 


5.440 


Fe 


26 


8.61 X 10-5 


0.058 


39.336 


9.914 


Ni 


28 


6.95 X 10-^ 


0.0054 


4.27 


10.662 








S = 1.58 


E = 61.2 





Table 2. Ionic species of iron at 20 MK 



Element 


Z-z 










K = Z'^ffX (keV) 


Fe XXll 


21 


21.98 


0.05 


0.43 X 10-' 


1.004 


6.57 


Fe XXlll 


22 


22.61 


0.14 


1.21 X IT^ 


3.152 


6.95 


FeXXIV 


23 


23.20 


0.25 


2.15 X 10-5 


6.232 


7.32 


FeXXV 


24 


2i.n 


0.56 


4.82 X 10-5 


15.381 


7.68 



with the m summation over m > Zixle)^^^, since recombination 
to level m yields only photons of e > Z^xlm^- If the source were 
so hot that all atoms were almost fully ionised the total for all Z 
would be, in this approximation, 

32;r r,V ^ ^ 1 F(e-ZW) 



3 Via 



for element abundances Az, with the same m summation limits. 

In reality even super-hot coronal flare temperatures are not 
high enough to equal the ultra-hot T » 10^ K needed to almost 
fully ionise all elements into their hydrogenic states, especially 
Fe, which is crucial in having by far the highest value of AzZ^ 
- see Table 1. Consequently, to deal accurately with jr for real 
flare data, we would have to take into account the actual ionisa- 
tion state of the flare plasma, which varies with time and location 
(being radically different in loop tops from loop footpoints), and 
actual forms of Qr{Z), Vz for non-hydrogenic ion stages. 

For our purpose of making first estimates we make the fol- 
lowing simpUfying approximations: 

- We treat all ions using hydrogenic Equations (9) - (11) but 
with suitably chosen Zg/f so that 



'z = ZgffX ; Q. 



32n 



RZ 



Z'- y2 



3 Via " nv'^^ 



(12) 



where Zg/y makes allowance for screening and other non- 
hydrogenic effects. While this will be a rough estimate for 
some ions, such approximations are often quite satisfac- 
tory for suitable Z^ff (e.g. Hahn and Krstic 1994, Erdas, 
Mezzorani and Quarati 1993). Here we adopt Z^ff such that 
hydrogenic Equation (12) gives the correct value of Qrz as 
given by exact calculations such as those of Arnaud and 
Raymond (1992) for Fe, which is the most important ion in 
our analysis. Typically, for an element of atomic number Z 
in an ionic state with z bound electrons left, Z^yy is between 
Z-z and Z - z -H 1 

Noting that Qr oc Xjnv' we include here only recombination 
to m = 1 (in the sense of the lowest empty level of the ion 
- hydrogenic with Z = Z^yj - not of the atom). Higher m 
contributions are weaker, being oc 1/m^ though extending 



to lower energies with edges at Z^^^xl^- These should be 
included in quantitative data fitting. 
- We focus on situations where the emitting region is near 
isothermal and either quite cool, so that only low Vz ele- 
ment recombination matters, or very hot so that high Vz el- 
ements (mainly Fe) are dominant. The former are typically 
loop chromospheric footpoints (thick target) and the latter 
very hot coronal loops which are either at the limb with their 
footpoints occulted, or are so dense as to be coronal thick 
targets (Veronig and Brown 2004). 

Under these conditions. Equation (11) becomes 



7«(e) = 



32;r 



2 2 



3 VI, 



a 



rip-Lz^ffZ^ffAz,, 



(13) 



where Az,^^ = Az^z,ff with <l>Ze// the fraction of atoms of ele- 
ment Z in ionic state Z^y/. 

Note that, since there is no integration over E here, if F(E) 
contains a sharp feature at an electron energy £«, such as a low or 
high E cut-off, this will appear in the recombination contribution 
to the photon spectrum j(e) as a series of sharp features at photon 
energies e(m, Z, £,) = Et+ Z^^^;^f/m^ ; m = 1, oo for every ion Z 
present. The same is true for broad features like smooth bumps 
or dips. This is in contrast with the bremsstrahlung contribution, 
in which such features are smoothed out by integration over E. 
Thus, even if Jr <*; Jb, it may have an important eff'ect in infer- 
ring F{E) from j{e) since this essentially involves differentiating 
jie) (Section 5). 



2.4. Element parameters and flare plasma ionisation 

The heavy element correction for bremsstrahlung, ^b, is almost 
independent of ionisation state (since the bremsstrahlung cross 
sections for atoms and ions of the same Z are essentially the 
same), being ~ 1-6 for solar abundances. On the other hand 
^RZeff - ^tff^z^ff depends on the number of empty ion levels 
available for recombination. The importance of fast electron re- 
combination radiation thus depends on the state of ionisation of 
the plasma in which the fast electrons are moving, which is pri- 
marily a function of plasma temperature T. 
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In Table 1 we show the values of Z, Z^Az - Cbz, 'Zf'Az = 
^Rz, Vz for various elements/ions whose = Z'^Az is large 
enough to be significant, if the element is sufficiently ionised. 
With (rtz ~ 40 for FeXXVI, Fe is by far the most important 
if conditions are such that it is highly ionised. The kT where 
maximum ionisation of an ion stage is reached is typically of 
the order OAZ^^j-x to ^^ffX- Table 2 we show more detailed 
values for several stages of ionisation of Fe (XXII-XXV, i.e. 21 + 
to 24+) with the appropriate Az^jj - Az^ff'^z^n for each of these 
Fe ionic states for the typical coronal flare case of T = 2 x 10^ 
K. These are taken from Amaud and Raymond (1992) as are 
the actual ionisation fractions we adopt later (Section 4) for the 
temperatures of the real flare we consider. 

The radiative recombination coefficients give Z^ff, which 
differ slightly from the Z values, as mentioned in Section 2.3. 
For the 2002 April 14 event, to which we return later, the peak 
flare temperature was 19.6 MK, ~ 5% of the iron appearing as 
Fe XXII (Fe2i+), ~ 14% in die Fe XXIII (^6^2+) state, ~ 25% 
appearing as Fe XXTV and ~ 56% as Fe XXV. The respective 
Z,ff values are 21.98, 22.61, 23.20 and 23.77. 

Broadly speaking in typical flare/micro-flare conditions we 
can consider the foUowing T regimes: 

- At r < W K ('cold') even H and oflier low Vz,„ ions are 
neutral so » for all Z. This would be typical of very 
dense cool chromospheric thick target footpoints relevant to 
deeply penetrating electrons. 

- For 10^ < r < 10*^ K ('cool') the predominant elements 
ionised are H, O, Mg, Si giving Szfj?z ~ 15. This is most 
relevant to upper chromospheric dense warm plasma reached 
by moderate energy thick target electrons. 

- At r > 10^ K ('hot') Fe is weU ionised up to about Fe XXV 
giving Sz^RZ ~ 50. This is relevant to the hot 'coronal' loop 
regime, hence either to (i) typical upper (SXR) flare loops of 
moderate density (thin target) whose HXR emission is seen 
in isolation either by HXR spectroscopic imaging or volume 
integrated but with the cool footpoints occulted because they 
are over the solar limb; or (ii) cases of coronal thick target 
loops (Veronig and Brown 2004) where the upper loop den- 
sity suffices to stop the fast electrons colUsionaUy. 

3. Local (thin target) HXR spectra of f-f and f-b for 
power-law F{E) with cut-off 

3. J. Bas\c expressions for JbJr 

To estimate how the fast electron recombination juie) compares 
with bremsstrahlung juie), we first consider the commonly stud- 
ied case of a power-law with a low energy cut-off 



E>Ec, 



(14) 



where Fc is the total electron flux at £ > Ec. Then, from 
Equations (4) and (14), we obtain for f-f emission 



jB(e) 



X 1 ; e <Er 



(15) 



jRz,ff{e) - yo- I) 33/2^^ — ~~ 



E,. 



X 0; 



e>E,+Zl^fX 



e<Er+Z: 



where 



(16) 



(17) 



So the total for all relevant Vz,ff is 
3.2. Ratio of Jr to js 

For this truncated power-law case, the ratio of f-b to f-f emissiv- 
ity is 



jR(e) 
jBie) 



-^7^Z^^>(e-£,)/;f) 



a25W5), 



e(keV) ^Z^ff>(e-E,)/x) & 



1 



-S-1 



(19) 



eff^- 

0.25'LzcffAz,ffZljjle(k&N). In pure ionised 



while for f-b emission from an ion of effective charge Z,,//, 



where each term in the summation is zero at e < + Z?, 

For e » Ec, ^ 0.25Sz.,,Az,,,z4 
H (Sz^jjz = 1) this is only 2.5% at 10 keV. This rather small 
value of *P must be the origin of the conventional wisdom that 
f-b can be ignored compared to f-f emission at HXR energies. 
However, this notion neglects several crucial facts: 

- At high coronal flare temperatures, where all elements are 
highly ionised, in plasmas of cosmic chemical abundances, 
heavy elements are the main contributors to the AzZ'^ sum. 
For the extreme ultra-hot case of near-total ionisation of all 
Z, and for modern solar coronal abundances the Sz factor 
is « 61.2, mainly due to Fe as discussed in Section 2.4 - see 
Tables 1 and 2. Note that Fe coronal abundance, for example, 
has been assumed to be 2.9 times photospheric Fe abundance 
(Feldman et al. 1992). Even higher factors of about 4 have 
been suggested (Dennis, personal communication). 

- At lower e the contribution from each Z^ff rises steeply to a 
sharp recombination edge aie = Ec + Vz, where the flux can 
be large, especially if Ec is small and 6 large. 

- At the edge, the [ ] factor in Equation (19) goes to [1 + 
^ffXlEcf^^ ■ This is because the flux of electrons emitting 
recombination photons of energy e is not the flux of those at 
£ > e, as for bremsstrahlung, but of those at £ = e - Z^^^jX- 
Consequently ^ is not negligible even at e » Ec- For fully 
ionised Fe alone, this factor is ~ 1 1 + lO/ii^ (keV)]*'^\ which, 
for 5 = 5 and at e = 10 keV, is 64, 11.4, 5.5 for Ec = 10, 
20, 30 keV respectively. Even for lower stage Fe ions (e.g. 
XXV), common in flare coronal loops, evidently recombina- 
tion must be a significant contributor to the HXR emission 
in those parts of the flare. 



3.3. Typical results in limiting regimes 

N.B. All spectrum figures in this paper (except Figure 5) have 
been plotted for a bin-width of 1 keV to match RHESSl's spec- 
tral resolution. However, in Figure 5 we use 0.01 keV resolution 
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Fig. 1. Actual shapes of non-thermal f-b and f-f spectra for different temperature regimes and non-thermal electron parameters. Note 
that the cool, hot and ultra-hot totals are almost identical and the dashed curves nearly indistinguishable for = 25 keV. 



so as to compare it with Plot A of Figure 4 to see how f-b edges 
would look if they were observed at a higher resolution. The 1 
keV binning smears out a lot of the edges of different elements 
that are clearly noticeable in Figure 5. Hence in Figures 3 and 4, 
the edges are not 'infinitely' steep as they should be; this is evi- 
dent in Figure 5 where they do look 'infinitely' steep due to the 
finer resolution. Also important to note is that the features seen 
in Figures 3, 4 and 5 are recombination edges and not spectral 
lines. None of the figures in this paper includes spectral lines - 
leaving them out shows more clearly where f-b edges exist in the 
HXR continuum. 

In Figure 1 we show for <5 = 3,5 the actual spectral shapes 
for Ec = 10,25 keV respectively in plasmas of normal solar 
coronal abundances, which are: ultra-hot (T » 10^ K; Fe is 
nearly fully ionised), hot (T = 2 x 10^ K; Fe well ionised up to 
Fe XXV) and cool (T - 10^ K; elements up to Si are almost fully 
ionised). In Figure 2 we show the ratios ^(e) for the ultra-hot, 
hot and 'cool' cases, respectively. The following key features of 
the hot thin target situation are apparent from these Figures: 

- The peak non-thermal f-b contribution, in each hot or ultra- 
hot case shown, adds at least 50% to the usual f-f one and in 
some cases (6 - 5,Ec = 10 keV) is up to 10 times greater 
(1000% increase) even when only ions up to Fe XXV are 
present. This is essentially due to the high abundance of Fe - 
much higher than thought when recombination spectra were 
first discussed (Culhane 1969, Culhane and Acton 1970). In 



appendix A we evaluate the efficiency with which f-b yields 
HXRs compared to f-f, and also derive the ratio T* for the 
case of a smooth F{E) with no cut-off. This proves, that in a 
hot enough plasma, far less electrons and power are needed 
than is found when only f-f is included and that, for smooth 
F{E), is largest for large 6 and low E spectral roll-over. 

- In the 'cool' case (T x 10^ K) of elements up to Si almost 
fully ionised, the f-b contribution is smaller but not in general 
negligible. For example, in the bottom left panel of Figure 1 
(5 = 5,£, = 10 keV), f-b is about 30 % of f-f at 15 keV 
energies. This is amply large enough to have a major impact 
on inferring F{E) by inversion or by forward fitting (Section 
5). 

- In hot plasma, Fe is by far the most important contributor of 
recombination radiation. 

- The peak ratio of f-b to f-f increases as 6 is increased and/or 
Ec is decreased. This is because f-b photons of energy e are 
emitted by electrons of energy E - V which have flux F{E - 
V) cc {E - Vy^ which is greatest when the minimum E - Ec 
is smallest, V is largest and the steepness S greatest. 

- Recombination edges are apparent for the elements with the 
highest values of Az^ffZ'^j:^ - Fe, Si, Mg and O and at energies 

e = Ec+Z^j:j:X, thereby creating the possibility of finding the 
location of a low energy cut-off Ec should one exist. 

- The harder asymptotic -y - 6 + I for f-f compared with y = 
6 + 2 for f-b (Equations (15) and (16)) results in an upward 
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Fig. 2. Photon flux ratio of non-thermal f-b to f-f emission for dififerent temperature regimes and parameters. Line styles have the 
same meaning as in Figure 1 . 



'knee' in the total spectrum clearly visible in Figure 1 for 
Ec = 10 keV but also present for higher Ec outside the e 
range of the Figure. This could be an important signature in 
data of a substantial f-b contribution. 



While the edge locations and the spectral shape trends will 
be roughly right, our use of the hydrogenic and Z^ff approx- 
imations, and adoption of unit Gaunt factors, mean that these 
curves/analytic forms can only be used for approximate quantita- 
tive fitting of real data. As far as we are aware (Kaastra, personal 
communication) the Gaunt factors, rates etc. have only ever been 
systematically evaluated for Maxwellian F(E) and sometimes 
for forms which can be written as sums of these (such as pure 
power-laws with no cut-off), and some occasional consideration 
of specific non-thermal spectra (e.g. Landini, Monsignori Fossi 
and Pallavicini 1973). Comparison of our Maxwellian results, 
in the unit Gaunt factor Kramers approximation, with those of 
Culhane for the same parameters shows the necessary correc- 
tions in the Maxwellian case to be significant for quantitative 
comparison with real data. In addition, in real cases the non- 
thermal emission will always be superposed on thermal contri- 
butions (especially important for the very hot plasmas of spe- 
cial interest here) and also in many cases on a thick target non- 
thermal contribution (unless this is from occulted footpoints), 
from the flare volume as a whole. In Appendix B we derive the 
generalisation of the above equations to the various cases in- 
volved in real flares, viz. finite volume thin targets, Maxwellian 



plasmas and thick targets for use in Section 4, where we evaluate 
the sum of all these contributions for a specific case. 

4. Some practical case study results derived from a 
real flare 

We saw above and in the appendices that the most favourable 
conditions for a substantial recombination contribution are when 
the maximum possible amount of the observable HXR source is 
a hot plasma (e.g. loop) at SXR temperatures. High density max- 
imises the emission measure but may make the source/loop colli- 
sionally thick and smear recombination edge spectral signatures 
of low energy cut offs. So an optimal case could be a loop which 
is just tenuous enough to be coUisionally thin and for which the 
cool dense thick target footpoints are occulted. (Footpoint re- 
moval by imaging is limited by RHESSI's dynamic range). Such 
sources will have a strong HXR source in the coronal loop. One 
such event was adopted as a basis for a case study, starting from 
the real event parameters. This was the 2002 April 14 event, 
which Veronig and Brown (2004) showed to be a hot, dense, col- 
lisionally thick loop with a strong coronal HXR source and no 
footpoints up to at least 60 keV. Thus the hot coronal source of 
non-thermal f-b emission was not diluted by cold footpoint thick 
target f-f emission though the f-b edges were smeared because 
the hot loop itself slowed the fast electrons to rest. In Figure 3 
we show the theoretical spectrum from a hypothetical resolved 
part of the coronal loop for two Ec values. We have evaluated 
the theoretical thermal, non-thermal and the whole volume hy- 
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Fig. 3. Spatially localised spectra from a loop with the 2002 April 14 event plasma parameters for two values of Ec. The left plot 
shows a very distinct iron edge at » 22 keV {- Ec + Vfe24+) and a much less predominant oxygen edge at « 15 keV (= Ec), whereas 
the second plot shows very distinct oxygen (ss 21 keV (= Ec)) and iron (« 28 keV) edges. This shows the value of recombination as 
an Ec diagnostic. The 'edges' appear to be of finite slope because of the finite (1 keV) resolution used. 



pothetical total Jsie), JR{e) (from Sections 2-3 and Appendix B) 
for such a loop, based on our approximate Kramers expressions, 
in three loop parameter regimes (Figure 4): 

- Plot A: With the actual hot thick target loop parameters 
found by Veronig and Brown, namely 6 - b.l; T - 19.6 
MK; L = 45 X 10*^ cm; A = 19.1 x lO'*- cm^; rip = lO" 
cm-3; = 4.9 X lO^" cm'^; = 5 X 10^^ sec"' above £i = 
25 keV. The total J is dominated by thermal f-b and f-f at 
low € but thick-target f-b at medium e and thick-target f-f 
at high e. Locally within the loop volume, if this were spa- 
tially resolved, the spectrum j would be like those in Figure 
3, where edges are clearly visible in positions corresponding 
to cut-ofF energies of 15 and 21 keV. At a higher resolution, 
these edges would look similar to the edges shown in Figure 
5. Should such edges be found in data, they can diagnose the 
all-important Ec parameter. 

- Plot B: With the actual parameters found by Veronig and 
Brown except with reduced by a factor of 25 so that the 
loop is collisionally thin above about 10 keV but with the 
footpoints hidden (limb occulted) so there is no cold thick 
target contribution. In this case the thermal emission is also 
much reduced because EM = Irt^AL is down by a factor of 
625. Somewhere between this and the first case should be the 
optimum condition for seeing maximum f-b contribution. 

- Plot C: The same as B but with the dominant cold footpoint 
thick target emission added to show its diluting effect. 

- Plot D: The same as C but with a reduced injection rate and 
so the thermal is more dominant than in C and this alters the 
total spectral shape a little bit. 

The upward 'knee' apparent in Figures 4 A,B at around 40 
keV due to the transition from a f-b to a f-f dominated spec- 
trum (cf. Section 3 and Figure 3) is rarely seen in data but may 
be present in some events (Conway et al. (2003)). A statistical 
survey of a large sample of events should shed light on con- 
ditions where non-thermal f-b is important. Also note that an 
upward 'knee' is present at the transition from a thermal- to a 



non-thermal-dominated spectrum. The position of this knee de- 
pends on the plasma temperature and may interfere with the f-b 
to f-f 'knee', which depends mainly on the Ec parameter. Hence, 
although for certain parametric conditions one may be able to 
notice two separate upward 'knees', if Ec is low and T is high, 
the 'knees' may occur at similar e and may not be distinguish- 
able in real data. 



5. The inverse problem - effect of f-f on F(E) 
inferred from data on j(e) 

We note again that, since even the thin target js involves an 
integral over E while Jr does not, any sharp features in F{E) 
would be smoothed out in the bremsstrahlung contribution to 
the photon spectrum but not in the recombination contribution. 
Consequently, an important way to study the effect of includ- 
ing f-b on the required properties of F(E) is to consider it as an 
inverse problem (Craig and Brown 1986) to infer F(E) from ob- 
served j(e). Here we consider the following experiment for the 
thin target case. (Thick target and thermal cases always involve 
even greater error magnification - Brown and EmsUe 1988). 
Generate the total 7(e) including f-b as well as f-f from a spec- 
ified Fi(E) and evaluate the F2(E) which would be erroneously 
inferred by solving the inverse problem ignoring the presence 
of the f-b term, as is currently done in all HXR data analysis, 
whether by inversion or forward fitting. 

By (4) and (11) the total f-f -l- f-b emission spectrum dJ/de 
from a homogeneous volume V can be written 

G(£)J£ + D2;^^^^^^z:^,Az.,,G(e-yz.,,), (20) 



where 



Hie)^ 



dJ 



Sar^ ^BfneC^ripV de 



e— ; G{E) = F{E)IE 



(21) 



and D is as given in Equation (A. 2). If we ignore the second 
(recombination) term in Equation (20), as has always been done 
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Fig. 4. The spectral components for 4 different hypothetical situations. We show these spectra by varying the parameters around the 
results in the Veronig and Brown 2004 paper that analyses the coronal thick target 2002 April 14 event. In all cases we keep the 
same values of 6 - 6.7, Ec = 10 keV and T = 19.6 MK. Plot A is for the thick-target coronal case with the actual event parameters 
Hp, Toe according to Veronig and Brown. Plot B was obtained for the same event parameters but with tip reduced 25 times to make 
the loop collisionally thin above 10 keV and with footpoint emission occulted. The injection rate is the same as Plot A so the density 
fraction of fast electrons is 25 times higher. The non-thermal emission is down by 25 times while the thermal is down by a factor 
of 625. Plot C is the same as B but with cold thick target footpoints included. The cold footpoint emission (motsly f-f) is dominant. 
Plot D is the same as C, but with an injection rate reduced by a factor of 25 so that the density fraction of fast electrons is the same 
as in Plot A. Evidently the detectability of the f-b contribution and of associated features in F(E) is sensitive to plasma parameters 
and observing conditions/geometry. 



in the past, for the Kramers f-f term, the inverse is just (Brown 
and Emslie 1988) 



G(e) = -H'{E). 



(22) 



The neglect of the second term can be thought of as an 'error' 
AH in our data and if we apply inversion formula (22) to this 
'data', ignoring the recombination 'error' we get a resulting error 
AG in the inferred G given by 



AG(£) = 



F2(E)-Fi{E) 



(23) 



■Z'„Az,,,G'(E-Vz,„). 



It is at once clear that any sharp change in j{e) i.e. in H(E), 
such as the presence of f-b edges, however small, can have a 
very large effect on the inferred F2{E). (If the inverse problem is 



addressed for more realistic smoother forms of f-f cross section 
than Kramers, the 'error magnification' is in general even larger 
- Brown and Emslie 1988, Plana et al. 2000). For a power law 
F with cut off around say 20 keV, analytically speaking this ex- 
pression gives infinite negatives in AG{E) at the spectral edges 
around 30 keV (for Fe). However when smoothed over a few 
keV and added to the f-f term the result would be a 'wiggle' 
in the F(E) solution in the 30-40 keV range. This is just where 
enigmatic features have been reported in some RHESSI spec- 
tra and variously attributed to the effects of photospheric albedo 
(Kontar et al. 2006), possibly pulse pile up (Plana et al. 2003), 
or a high value of Ec (Zhang and Huang 2004). 



Another case providing insight is that of a smooth shifted 
power-law G(E) - A(E + E,y^^\ which has no edges though 
the corresponding F(E) has a smooth peak at E = E,/6. In this 
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Fig. 6. Fractional error (AG/G) in G (Equation (25)) as discussed in Section 5 for £, = 5, 10,20 keV respectively for a shifted 
power-law due to inference of G from H ignoring the presence of recombination. 



case the fractional error in G due to applying (22) ignoring the 
recombination term can be expressed as 



AG(E) 
G(E) 



= (6+1) 



D 



E + E, 



1 



^^+2 



/(£ + £.) 



where each term in the Zpyy sum is zero for E < Vz^yy 

In the case of recombination onto Fe XXV alone (hot 
plasma), this gives for <5 = 5, 



AG 



lOkeV 
E + E, 



[1 -1 keV/{E + E,)]- 



(25) 



which is shown in Figure 6 for - 5, 10, 20 keV. Evidently 
errors due to neglect of recombination can be large at low E. 
The reason is that the Ze/f recombination contribution to the 
bremsstrahlung solution for G(E) at E comes from the slope of 



G, and not just G itself and at E - Vz^ff not at E. Figure 6 is 
similar to Figure A. 2 because F2/F1 = G2/G1 = 1 + AG/Gi. 
This error has very serious consequences for past analyses 
/^A-pf HXR flare spectra, at least in cases where a significant hot 
' dense coronal loop is involved. For example, the f-b emission 
spectrum is most important at lower energies (5-30 keV or so), 
depending on the plasma temperature T and low energy elec- 
tron cut-off or roll-over Ec, £* and is steeper than the free-free. 
This will off'set some of the spectral flattening caused around 
such energies by photospheric albedo (Alexander and Brown 
2003, Kontar et al. 2005) resulting in underestimation of the 
albedo contribution and hence of the downward beaming of the 
fast electrons. This fact would weaken the finding of Kontar and 
Brown (2006) that the electrons are near isotropic, in contradic- 
tion of the usual thick target description, but for the fact that the 
flares they used had rather hard spectra and substantial footpoint 
emission - conditions where the f-b correction should be rather 
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small. Nevertheless it illustrates that care is needed to ensure f-b 
emission is properly considered. 

Finally, recognising the presence of the f-b contribution, 

one can in fact convert integral Equation (20) into a differen- 
tial/functional equation for F{E) by diflFerentiation, namely 

G(E) - i)2z.,,>£/;,"^^z,,,z4^^G'(£ - Z^^ffX) = -H'(E), (26) 

which is a whoUy new class of functional equation in need of 
exploration. 

6. Discussion and Conclusions 

It is clear from our findings that ignoring non-thermal f-b contri- 
bution as negligible, as has been done in the past, is erroneous. 
Even if we ignore coronal enhancement of element abundances, 
and use photospheric abundances, f-b contribution can be very 
significant. In certain flaring regions, especially in dense-hot 
coronal sources or occulted loop-top events, fast electron recom- 
bination can be of vital importance in analysing data properly 
and in inferring electron spectra and energy budgets. It can have 
a major influence on inferred electron spectra both as an inverse 
problem and also in forward fitting parameters, including the im- 
portant potential to find and evaluate low-energy electron cut- 
ofi^s, which are vital to flare energy budgets. While incorporat- 
ing f-b into spectral fitting procedures wiU make it considerably 
more complicated, an advantage is that the f-b, unlike the f-f, 
contribution retains its J{e) signatures of any sharp features in 
F(E). 

A major consequence of the low energy f-b contribution is 
that, to fit an actual photon spectrum, less electrons are needed, 
than in f-f only modelling, at the low E end, which is where 
most of the power in F(E) lies. For example, if we consider the 
case 6 = 5,Ec = 10 keV and ionisation up to Fe XXV, then we 
see from Figures 1 and 2 that inclusion of f-b increases j by a 
factor of 2-10 in the 15-20 keV range for 6 = 3-5. Thus, to get 
a prescribed j in that range we need only 10 - 50% as many 
electrons as inferred from f-f emission only. 

We also note that the importance of non-thermal f-b emission 
is greatest when non-thermal electrons are present at low E and 
with large 6 such as in microflares with 'hard' XRs in the few 
to ten KeV range (Krucker et al. 2002). Such low energy elec- 
trons have short collisional mfps and so are more likely to emit 
mainly in hot coronal regions, if accelerated there. Microflares 
are therefore important cases for inclusion of f-b. 

Before we conduct any precise fitting of F(E), involving the 
f-b contribution, to real data (e.g. from RHESSl) and include it 
in software packages it wiU be important to include, for both f-b 
and f-f, more accurate cross-sections with Gaunt factors etc. and 
ionisation fractions as functions of plasma temperature. By do- 
ing this, it will be possible to show, for certain events, how vital 
recombination is and to improve our understanding of electron 
spectra and their roles in flares. However, our Kramers results 
already bring out the fact that recombination should not be ig- 
nored in the future, and that it may be invaluable in some cases 
as a diagnostic of the presence or otherwise of electron spectral 
features. 
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Appendix A: Efficiency and smooth F(E) 

A.1. Comparison of the efficiency of f-b versus f-f HXR yield 

In Section 3.1 and 3.2, we predicted the j'b, from a power- 
law F(E) and found that the Jr contribution could sometimes be 
more important than j'b. It is of interest therefore, to consider 
the following question. If one observes a power-law /(e) oc e"'' 
above some e > Vz, what electron flux FniE) would be needed 
to generate it in a plasma of solar abundances purely by non- 
thermal electron recombination on ion Z<,// as compared with 
the Fb{E) required to do so purely by f-f bremsstrahlung? If we 
write, from Equation (15), the latter as Fb(E) = CE~^^^ then the 
former has to be, by Equation (16), 



Fr{E) = C{E + Vz^^^r^/DZ^^ffAz^^ 



where 



D = 



0.04keV 



(A.1) 



(A.2) 



and the ratio measuring recombination efficiency relative to 
bremsstrahlung is 



FsiE) 
Fr{E) 



= yZ%fAz,„^\l + Vz,„IEy, 



(A.3) 



which we show in Figure A.l for y = 5 in terms of each of 
the dominant f-b contributions from fully ionised O, Mg, Si and 
Fe respectively while the f-f is for all elements. Evidently non- 
thermal recombination could be dominant over bremsstrahlung 
up to many 10s of keV as the most efficient HXR source if 
the electrons are emitted entirely in a plasma hot enough (7" « 
20MK) for elements up to Fe 24-1- to be ionised and is significant 
even at lower temperatures. 

In terms of the total required electron fluxes Fm, Fb\ above 
energy £i, the ratio is 



Fbi 
Fri 



— = ^Z\Az„T-i^ + Vz„/EiV-^ 

J7„. y-2 eff '^'itE,'- ^'ff' 'J 



0.02Z4 4z, 



10 keV 



[1 + Vv„IEi\ 



y-1 



(A.4) 



which is about 10 for Fe, 0.25 for Si and 0.1 for Mg and O at 
£i = lOkeV. 

At higher electron energies (E >v 17 keV), O becomes more 
efficient than Mg, as can be seen in Figure A.l, because of the 
combined effects of the AzZ* factor and the term containing Vz- 

A.2. Ratio of jR to jb for an example of a smooth F{E) with 
no cut-off 

All of the above results are for F{E) with a sharp cut off Ec. To 
illustrate how the appearance of j(e) is modified by inclusion of 
f-b as well as f-f for a smooth F{E), a simple case to evaluate is 
F{E) oc E{E + E,)-^-^, which behaves as E'^ at £ » £. but has 
a smooth roll-over at E^jS. It is simple to show that the resulting 
jsie) cc {E + EtT^/S for f-f alone and that the ratio of f-b to f-f 
in this case is, for ion Ze/f alone. 



^smooth ~ 



e + E, 



1- 



e + Ef 



-6-1 



(A.5) 



which is shown in Figure A.2 for 5 = 5, Zg/f = 23.77 and = 
5, 10, 20 keV. We see again that ^smooth is largest for large 6 and 
for small 
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Fig. A.2. The ^smooth as discussed in Appendix A and Equation (A. 5). It is the ratio of Jr to Jb for the smooth F(E) oc E(E + £») 
for £, = 5, 10, 20 keV respectively. 
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Appendix B: Whole Flare Thin Target, Thermal, and 
Thick Target Expressions for f-f and f-b HXR 
Emission Spectra 

Here we extend the above results on local emissivities j{e) to 
estimate total spectral emission rate J(e) (photons sec"' per unit 
e) from extended flare volumes as required for real flare data. 



ii] 

1; 



and by Equation (16), 



-s 



e < E, 



(B.l) 



B.1. Thin Target Coronal Loop 

A thin target is one in which F{E) is not significantly modified 
by energy losses or gains over the volume. For a loop of half 
length L, transverse area A, volume 2AL and density Wp, the total 
emission rate spectra contributions Jthini^) are for a power law 
F{E) with a low energy cut-off, by Equation (15), 



7«„,„(e)= (5-1)^^. 



]X 2n„ALF, 



X 0; 



El 

-5-\ 



Xl.7 



e>E,+Z: 



(B.2) 



Jb„,,„{^) — ^ «i,c inpALF, ^ scalcd by the plasma volume. 



where the summation is over all Z^jj < (e-Ecl These spectral 
shapes 7(e) are of course just the same as the thin target j forms. 
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Fig. B.l. Non-thermal f-f and f-b spectra for the thick target case (Equations (B.l 1) and (B.12)) shown for 2 different temperatures: 
20 MK that is pertinent to events such as the 2002 April 14 event and 10 MK, which is more in the range of 'microflare' temperatures. 
It is interesting to note the three distinct energy regimes for the f-b spectrum, namely: e < Vfe', Vfe < e < Vfe + Ec', e > Vfe + Ec- 
Clearly f-b is very important in the 10-50 keV range, precisely where albedo issues are also important. 



B.2. Hot Coronal Loop Thermal Emission (In the Kramers 
approximation) 

Both f-f and f-b emissions are included in the standard analy- 
ses (e.g. Mewe et al. 1987, Dere et al. 1996) of isothermal hot 
plasma contributions to flare spectra, using full cross sections 
and ionisation balance expressions. It is therefore surprising that 
f-b is omitted from calculations of non-thermal emission, espe- 
cially at low e, where electrons of comparable energy are present 
in both thermal and non-thermal populations. In applying our 
study of the non-thermal f-b to real data we wish to include ther- 
mal emission as it is important at energies under about 20 keV 
and so dilutes the visibility of non-thermal contributions. In or- 
der to treat the thermal and non-thermal j consistently and allow 
meaningful comparisons we use the expressions for the thermal /' 
relevant to the Kramers cross sections just as in the non-thermal 
case - but see remarks previously and below concerning Gaunt 
factors and absolute accuracy of our results. 

For an isothermal plasma the local Maxwellian electron flux 
spectrum is 



fherm' 



(E) 



nnie 



1/2 



rip exp(-E /kT), 



(B.3) 



which, by Equation (4), gives for the thermal bremsstrahlung 
emission from a uniform loop 



J 



B therm 



(e) 



I6an 



1/2 InlALe-'l'''^ 



e(kTy/^ 



(B.4) 



and for the recombination 



jRtherm(e) — 



exp 



2n-„AL 



(B.5) 



These results can be compared with those of CuUiane (1969) 
and Culhane and Acton (1970) who were among the first to ex- 
plicitly address the X-Ray spectrum from hot coronal plasmas. 



Using the Kramers cross sections is essentially equivalent to set- 
ting to unity all Gaunt factors in their expressions. When we 
do so, the £,T dependences of our Jutheim, JBrherm are identical 
to theirs - e.g. JRthennI JBthenn IS independent of e, the only differ- 
ence being that our JR,herm is much larger (in absolute value) than 
theirs, mainly because they used the very much lower value of 
Az for Fe believed at that time. Examination of the e, T depen- 
dences of Culhane's Gaunt factors shows that they affect quite 
significantly both the f-f and the f-b spectra from a Maxwellian 
F(E) and we should expect the same to be true for non-thermal 
F{E) like power-laws. Thus, any accurate absolute comparison 
of predictions with data will require incorporation of appropriate 
g, G. However, these do not affect the absolute orders of magni- 
tude of JRthenn, JBthenn nor the dependencies on Up, V, Fc etc., nor 
the locations of edges. So, for the present purpose of demonstrat- 
ing the importance of f-b, the Kramers expressions will suffice. 

B.3. Thick target (dense loop or footpoint) f-f and f-b 
emission spectra 

In the thick target case, j evolves in space along with the energy 
losses of the electrons. To find j locally one uses the continuity 
equation (Brown 1972) and then integrates over volume to get 
J. However, to get the whole volume J, it is actually simpler 
(Brown 1971) to start with the electron injection rate spectrum 
ToiEo) electrons/sec per unit injection energy Eg and use the 
expression 



Jihickie) 



)r]{e, Eo)dEo, 



(B.6) 



where //(e, E,,) is the total number of photons per unit e emitted 
by an electron of energy Eg as it decays in energy. For purely 
collisional losses dE/dN - -K/E with K - Ine'^A, e being the 
electronic charge and A the Coulomb Logarithm. Then 



Eo) 



— f E—dE 
KJe de 



(B.7) 



for the relevant radiation cross section dQjde. Note that this as- 
sumes H to be uniformly and fully ionised along the electron 
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path. For partially ionised H the energy loss constant K is re- 
duced but this situation is not relevant to our hot source situa- 
tions. 

For our Kramers dQ/de f-f and f-b expressions (3), (7) and 
(9), the resulting expressions, in the case where Az^^^ are uniform 
along the path, Equation (B.7) gives 



0; e>Eo 



(B.8) 



and 



T]Rzie,Eo) = 



1; 
0; 



Ke 



Eo>e + Z'^ffX 

Eo<e + Zlff. (B.9) 
For a power-law injection rate spectrum of spectral index 6o, 



VIZ 



ToiEo) = i6o - 1)|^ 



(B.IO) 



where Toe is the total rate above low energy cut-off £oe> the 
expressions for the non-thermal emission spectra are then by 
Equation (B.6) 



'Bthick 



(e) 



3 (6„-\)(S„-2)K 
\-<5»+l 



[(^«-l)|-(^»-2)]; e<E, 



(B.ll) 



and, for ion Z^//, 



jRZ,ffthicki^) — 



"tC r X Toe 
33/2^ WZefi Ke E„c 

■<5<,+l 



e>E„r+Z: 



AffX <e<Eoc + 



0; 



e<Z; 



(B.12) 



For the case of a cold thick target footpoint the total can be 
almost as small as 1 if only hydrogen and some low elements 
are ionised and even zero if T < 8000 K or so (there being almost 
no charged ions present). In these sources the f-b contribution is 
negligible or at most a very small correction. For a colUsonally 
thick hot loop ^r is, however, very much higher. 

The main distinction of these hot thick target spectra com- 
pared to hot thin targets is that the decay of all electrons to zero 
energy means that the signature of the cut off Egc in the injection 
spectrum appears not as a discontinuity in /(e) but only in its 
gradient J'{e). This gradient break is very noticeable in Figure 
B.l at energy e = Ec + Vpe- So, even in the thick target case, 
spectral diagnosis of any Egc present is possible. The recombina- 
tion edges themselves appear at the relevant ionisation energies 



e = Veff, these being from thick target electrons decelerated to 
zero E. These non-thermal recombination spectral edges are then 
down in the energy regime below 10 keV which is complicated 
by Fe lines etc., making the interpretation of To there, and of the 
lines, more difficult. 
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